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Abstract. We study the systematic problems in determining the local dark matter den- 
sity Pdm(^o) from kinematics of stars in the Solar Neighbourhood, using a simulated Milky 
Way-like galaxy. We introduce a new unbiased method for recovering Pd,„(i?o) based 
on the moments of the Jeans equations, combined with a Monte Carlo Markov Chain 
(MCMC) technique and apply it to real data [ 1 j. 



1 Introduction 

A good constraint on the local dark matter (DM) density gives us useful information about local 
Galaxy dynamics. It is also vital for direct DM search experiments. With near-future surveys promising 
a dramatic improvement in the number and precision of astrometric, photometric and spectroscopic 
local stars data, it is timely to revisit the systematic problems arising in determining the local matter 
density, as well as the local dark matter density - pdm(^G) - from the vertical motion of stars at the Solar 
Neighbourhood. To do this, we use a high resolution N-body simulation of a Milky Way-like galaxy 
to test the standard methods in the literature. We also introduce a new method for recovering pdm(^o) 
based on the moments of the Jeans equations, combined with a MCMC technique to marginalise over 
the unknown parameters. This method relies on very few assumptions, but requires more precise data 
than the classical methods adopted in the literature. Given sufficiently good data, our method recovers 
the correct local dark matter density even in the face of disc inhomogeneities, non-isothermal tracers 
and a non-separable distribution function 1 1 1. 

2 Method and results 

The standard methods to determine the local density from the kinematics of stars II2I3I4I5II rely on 
several assumptions: 

1. The disc is in equilibrium (steady state assumption) M2I3I4I5I ; 

2. The local dark matter density is a constant pdm(^G) in the z-range considered (<1 kpc) II2I3I5II : 

3. Tilt terms involving cr^s in the Jeans equation are negligible |2 3 4 5|; 

4. The local visible matter is a sum of isothermal components [1213 15 1 ; 

5. The rotation curve of the Milky Way (MW) is flat li2i3ii4i5J : 

6. The distribution function of the tracers is a function only of the vertical energy: / = /(£/) 1213151 . 
In this case the density law of a stellar tracer v(z) can be predicted from its vertical velocity distii- 
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' Actually, this implies that a^R = 0, rhather than just very small. 
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Fig. 1. Density of models explored by the MCMC (for 8 volumes around the disc) represented as shaded areas of 
different colours (blue=pdn,, red=ps). The filled dots represent the corresponding actual values from the simulation 
(with Poisson error bars), the horizontal segments are the 90% confindence ranges. Left panel: standard method, 
Right Panel: MA method. 



The comparison between the predicted and the observed v(z) of the tracer stars constrains the local 
potential, and thus the local visible and dark matter densities f2'3'51. Our new method, based on the 
solution of the Jeans-Poisson system, (i) relies on a Minimal set of Assumptions (MA method [ 1)) - 
namely[r||2][3]- and (ii) uses a MCMC technique to marginalise over unknown parameters. The former 
makes the method - given good enough data - robust to model systematics; the latter allows us to cope 
with incomplete or noisy data and model degeneracies. In the MA method the density law of a tracer 
population v{z) - in a given potential - is predicted from the solution of the Jeans equation 

v(0) <r?(:) '^^P 1^ Jo a] dz "^^) 

where <t^{z) is the vertical velocity dispersion law. 

To test both the standard and the MA method, we build a high resolution equilibrium N-body 
model approximating the Milky Way that satisfies all of the usual assumptions. We then evolve the 
disc over ~ 4 Gyr to form a complex disc structure with a strong bar and spiral waves similar to the 
MW. We fit v{z) to find the best value of of the local visible and dark matter densities ps and pdm at 8 
different angular positions around the disc, representing the Solar Neighbourhood (at /? = 8.5 kpc from 
the Galactic centre). The application of the standard and the MA methods to the unevolved simulation 
gives similar results and recovers the correct value of ps and pdm- However, we find that we need data 
extending up to ~ 3 times the half mass scale height of the disc to break a very strong degeneracy 
between ps and p^^ (see also |5|). The situation is very different for the evolved disc that does not 
satisfy assumptions|4]|5]and|6] While it is easy to correct for|4]and|5] the failure of hypothesis |6]makes 
the standard methods fail in recovering the coiTect ps and pdn, (see fig.[l]i. 

We illustrate our MA method by applying it to Hipparcos data |2 3|, extending to ~ 700 pc. We 
first make the usual assumption that the (A, F and K giant) star tracer populations are isothermal. This 
recovers p^m - 0.003;^Q Q[j^ Mopc"-' (with 90% confidence), consistent with previous determinations. 
However, the vertical dispersion profile of these tracers is poorly known. If we assume instead a non- 
isothermal profile similar to the blue disc stars from SDSS DR-7 |6|, we obtain an equally good fit to 
the tracer density profile, but with ut with pd^ = 0.033^g QQg M© pc^^*. This highlights the fact that it is 
vital to measure the vertical dispersion profile of the tracers to recover an unbiased estimate of pdm. 
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